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The femtosecond optical pump-probe technique was used to study dynamics of photoexcited 
electrons and coherent optical phonons in transition metals Zn and Cd as a function of temperature 
and excitation level. The optical response in time domain is well fitted by linear combination of 
a damped harmonic oscillation because of excitation of coherent E2 g phonon and a subpicosecond 
transient response due to electron-phonon thermalization. The electron-phonon thermalization time 
monotonically increases with temperature, consistent with the thermomodulation scenario, where at 
high temperatures the system can be well explained by the two-temperature model, while below w 
50 K the nonthermal electron model needs to be applied. As the lattice temperature increases, the 
damping of the coherent E^ a phonon increases, while the amplitudes of both fast electronic response 
and the coherent E^ a phonon decrease. The temperature dependence of the damping of the E^g 
phonon indicates that population decay of the coherent optical phonon due to anharmonic phonon- 
phonon coupling dominates the decay process. We present a model that accounts for the observed 
temperature dependence of the amplitude assuming the photoinduced absorption mechanism, where 
the signal amplitude is proportional to the photoinduced change in the quasiparticle density. The 
result that the amplitude of the E^ g phonon follows the temperature dependence of the amplitude of 
the fast electronic transient indicates that under the resonant condition both electronic and phononic 
responses are proportional to the change in the dielectric function. 



PACS numbers: 78.47.+p, 78.20.-e, 63.20.Kr, 78.66.Bz 

I. INTRODUCTION 

Dynamics of nonequilibrium electrons and phonons in 
metals, semiconductors, and superconductors have been 
the focus of much attention because of their fundamen- 
tal interest in solid-state physics. In particular, electron- 
phonon interaction is the center of the interest for under- 
standing the nature of nonequilibrium electron relaxation 
[electron-phonon (e—ph) thermalization], of dressed elec- 
trons (an electron plus its phonon cloud) , of phonon self- 
energy (real part corresponds to the change in phonon 
frequency), and of phonon-induced electron-electron in- 
teraction (this effect gives the birth of the Cooper pair 
in superconductors). [lj 

In metals, relaxation dynamics of optically excited 
nonequilibrium electrons has been extensively studied 
by transient-reflectivity (TR) or transient- transmis- 
sion (TT) pump-probe techniques using femtosecond 
lasers. 0,03, EL III 11 The e—ph thermalization occurs on 
a subpicosecond time scale in metals and has been tra- 
ditionally described by the use of the two temperature 
model (TTM).|E[1|1I3 Besides the observation of the 
dynamics of nonequilibrium electrons, coherent acoustic 
phonon pulses have been investigated by using the pump- 
probe angular deflection technique. |Tot ITU The acoustic 
phonon pulse is generated through transient heating of 
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the sample surface induced by excitation with intense 
laser pulses. The response from acoustic phonon pulses 
is typically observed on a several hundred picosecond 
time scale (or gigahertz frequency range) [l^. Since the 
electron-phonon thermalization that we are considering 
takes place on a subpicosecond time scale (or terahertz 
frequency range), Q coherent optical phonons in metals, 
which arc impulsively excited and decay within a few pi- 
coseconds, may play an important role in relaxation of 
hot electrons. 

Recently, Melnikov et al. have observed the surface 
coherent optical phonon (u = 2.9 THz) at the surface of 
Gd metal using the second harmonic generation (SHG) 
technique. 13] They focused their attention into the gen- 
eration mechanism of the coherent optical phonon and 
coupled coherent magnon. The surface coherent optical 
phonon was excited by a transient charge separation at 
the surface, which is similar proposal to the displacive 
excitation mechanism (DECP).[lJ] Furthermore, Boven- 
siepen et al. studied both surface and bulk coherent 
optical phonons in Gd metal using the SHG and TR tech- 
niques, respectively. They discussed the transient fre- 
quency shifts observed for the both phonons. Watanabe 
et al. revealed coherent surface vibration on Cs/Pt(lll) 
system by using the SHG technique. 16] The coherent 
stretching vibration (to = 2.2 THz) was generated by res- 
onant impulsive stimulated Raman scattering (ISRS).[T^] 
and the damping of the coherent surface vibration was 
dominated by pure dephasing caused by hot electrons at 
absorbate. Thus, the SHG technique has mainly enabled 
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studying surface coherent optical phonons on metallic 
surfaces. Nevertheless, investigations of bulk coherent 
optical phonons in metals using conventional TR or TT 
pump-probe techniques are still few mainly because of 
the very short optical penetration depth in metals stem- 
ming from the absorption by free electrons (Drude ab- 
sorption). Therefore the change in the refractive index 
associated with the coherent phonon oscillations are very 
weak in the bulk, making it very difficult to measure. Al- 
though study of the surface coherent optical phonon in 
metal surfaces offers interesting applications, e.g., con- 
trolling chemical reactions on metal surfaces, there is 
also tremendous motivation to observe and control the 
bulk coherent optical phonon in metals in order to study 
transient melting as a precursor of phase transition[l£(, 
as well as controlling electron-phonon scattering. [19j 

In absorbing media, e.g., in semimetals, resonant ISRS 
has been proposed as the main generation mechanism 
of the coherent Ai s phonon. [l7l Ecj] Stevens et al. con- 
sidered both the standard Raman susceptibility and the 
electrostrictive tensor as the driving force for resonant 
ISRS. 20] Since electrostrictive force is proportional to 
the dielectric function, the amplitude of the coherent 
phonon depends on the incident photon energy. |20j In 
some metals, interband transition occurs in the near in- 
frared region, which appears as a peak in the imagi- 
nary part of the dielectric function, and the interband 
transition could enhance the amplitude of the coherent 
phonons through the resonant ISRS process. Therefore 
the realization of generation and observation of bulk co- 
herent optical phonons in metals using the TR pump- 
probe technique requires (i) existence of k = optical 
phonon, which can be observed by Raman scattering and 
(ii) a spectral peak in the imaginary part of the dielectric 
function, which is related to interband transition. Appro- 
priate candidates of the metallic samples that meet the 
above requirements, are Zn, Cd, and Mg, in which k = 
optical phonon modes were, in fact, observed in Raman 
scattering. [21) 

In this paper, we report on the observation of both bulk 
coherent optical phonon and nonequilibrium electrons in 
Zn and Cd using a femtosecond pump-probe reflectivity 
technique with high sensitivity of AR/R = 10~ 7 . The 
coherent Ei g phonon, whose oscillatory frequency is in 
the terahertz region, is generated by optical excitation 
of nonequilibrium electrons (interband transition). The 
amplitude of the oscillations decreases dramatically as 
the temperature is raised from 7 K to room temperature. 
Importantly, the temperature dependence of the ampli- 
tude of oscillation closely follows the temperature depen- 
dence of the subpicosecond electronic transient, suggest- 
ing that under the resonant condition both electronic 
and phononic responses are proportional to the change 
in the dielectric function. The relaxation time of the 
fast electronic response decreases with decreasing tem- 
perature, showing saturation at temperatures below 100 
K. In fact, this behavior closely follows the temperature 
dependence of the electron-phonon thermalization time 
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FIG. 1: (color online) Schematics of femtosecond pump-probe 
experiment. BS: beam splitter, A/2: half wave-plate. The 
difference of the photocurrent of the two Si-PIN detectors 
was amplified, and thus reduced current noise from the laser 
source. 



observed in simple metals, such as Au, Ag |9j]. In order 
to explain the peculiar temperature dependence of the 
amplitude of the electronic transient, we have derived 
a simple model, where the photoinduced change in re- 
flectivity is governed by photoinduced absorption, where 
the initial states lie near the Fermi energy. Assuming the 
so-called two-temperature model and energy conser- 
vation law we were able to reproduce the temperature 
dependence of the transient amplitude extremely well - 
further supporting the association of a fast transient with 
electron-phonon thermalization and correlation between 
electronic and phononic responses via the dielectric func- 
tion. We have extended the range of pump fluences up to 
6 mJ/cm 2 and study phonon softening as a spectator of 
thermal expansion. In Sec. V, we discuss the tempera- 
ture dependence of the amplitude of the coherent optical 
phonons in Zn in terms of Stevens' model. |20j 



II. EXPERIMENTAL TECHNIQUE 

The samples used here were single crystals of Zn and 
Cd with cut and polished (0001) surfaces. The femtosec- 
ond pump-probe measurements (Fig. 1) were carried out 
in a temperature range from 7 to 295 K using a cryostat. 
The light source used was a mode-locked Ti:sapphire 
laser whose pulse width was ~ 20 fs as determined by 
autocorrelation of ~ 30 fs,|2^ and the repetition rate 
of 87 MHz. The pump-beam power was varied from 20 
to 120 mW, while the probe-beam power was fixed at 5 
mW. The polarizations of the pump and probe beams are 
orthogonal to each other to avoid picking the scattered 
pump beam on the detector. We used a lock-in detec- 
tion where the pump-beam was mechanically chopped at 
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2 kHz. Both pump- and probe-beams were focused onto 
a diameter of ^ 100 /im on the sample. The penetra- 
tion depth of the laser light with a wavelength of 800 nm 
(= 1.55 eV) was estimated to be ~ 13 nm based on the 
absorption coefficient of 7.59xl0 5 cm -1 . Therefore the 
contribution to the signal from the surface oxide layer is 
negligibly small. At the pump fluence of 9.2/zJ/cm 2 (60 
mW) we estimated the maximum electron temperature 
rise AT™ (assuming a thermalized electron gas) to be < 
100 K at 7 K.[23| Thus the experiments can be considered 
as being in the reasonably weak excitation regime. H,0| 
the TR signal (AR/R) was measured to extract the relax- 
ation dynamics of noncquilibrium carriers and the bulk 
coherent optical phonon as a function of the time delay 
after excitation pulse, pfij 



III. GENERATION OF COHERENT OPTICAL 
PHONONS 

Traditionally, generation of coherent optical phonons 
has been described by either DECP^4| (for absorbing 
media) or ISRS 26] (for transparent media). In DECP 
mechanism, photoexcitation induces changes in the elec- 
tronic energy distribution function, and as a result, the 
crystal lattice begins to oscillate around the new equilib- 
rium position; Qo(t) — nn(t), where k is a constant and 
n(t) is the photoexcited carrier density. [T^j In the first or- 
der only the Ai fl totally symmetric modes are coherently 
excited by the DECP mechanism. Recently, Stevens et 
al. have suggested that the combination of DECP and 
ISRS is possible in semimetals. We briefly introduce their 
model below. |20j 

The impulsive stimulated Raman scattering equation 
of motion for an optical vibrational mode whose normal 
coordinate is labeled Q is given by[2(], H(j 



force F(t) is expressed as follows :[2C 



d 2 Q dQ o 



k,l 



where T and Qo are the vibrational damping and the 
frequency, respectively, N is the number of cells, v c is 
the volume of the unit cell, while E^i are components 
of the electric field delivered by the broad spectrum of 
the ultrashort pump light. According to wave- vector 
and energy-conservation rules, a phonon mode with fre- 
quency u>i — uJk = and wave- vector ki — kk = fco can be 
excited. [2a, l2o] The Raman tensor Rki denotes both the 
standard Raman tensor Xm an d the electrostrictive force 
7Tjy. Equation (1) is valid only if Rj~i does not depend on 
frequency. In absorbing media, such as semimetals and 
metals, the Raman tensor Rki depends on frequency and 
F(t) becomes a function of the frequency. According to 
the theoretical work for the electrostrictive force irj^ as 
the generation mechanism, assuming that the dielectric 
function s(u>) varies slowly within the spectral width of 
the pump pulse, the Fourier component of the driving 
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e lilt \E{t)\ 2 dt. (2) 



Here E(t) is the electric field of the pump pulse, i.e., 
|-E(i)| 2 corresponds to the intensity of the pump light 
\E(t)\ 2 = I(t). The real part of e leads to F(t) oc \E(t)\ 2 , 
which is impulsive and gives Q(t)oc sinf2ot, whereas the 
imaginary part of e leads to F(t) oc J_ \E(t')\ 2 dt' , 
which is displacive in character and gives Q(t)oc cosfiot. 
Since |dRe(e) / duj\-^Jm(e) / TIq with fio being the phonon 
frequency, imaginary part of the dielectric function domi- 
nates the driving force. Thus we expect that the coherent 
optical phonon generated through resonant ISRS exhibits 
a displacive behavior. Solving Eq. (1) for the undamped 
harmonic condition (T = 0) yields the coherent phonon 
amplitude. p0| 



Im(e)Nv c ^ 

4nhn 2 



e lUot \E(t)\ 2 dt, 



(3) 



where So is the deformation potential. According to Eq. 
(3), the amplitude of the coherent phonon depends on 
the imaginary part of the dielectric function as well as 
the pump intensity I(t) if the temperature is constant, so 
that thermal expansion and anharmonicity of the crystal 
lattice, which result in modification of v c and f2o, are 
negligible. 



IV. EXPERIMENTAL RESULTS AND 
ANALYSIS 

Figure 2 shows the time- resolved TR signal (AR/R) 
obtained for Zn and Cd at 7 K with the pump fluence 
of F pump = 9.2/iJ/cm 2 . The response of the reflectivity 
change consist of two components. One is the initial tran- 
sient nonoscillatory response because of excitation and 
relaxation of nonequilibrium electrons (fast electronic re- 
sponse), which decays in a few hundred femtoseconds. 
Since the interband electronic transition near the L-point 
occurs at around 800 nm (1.55 eV),|23 this interband 
transition dominates generation of nonequilibrium elec- 
tron distribution in Zn. 

The second component is the oscillatory signal due 
to generation of the coherent lattice vibration (coherent 
phonon response) . The time period of the observed co- 
herent oscillation is several hundred femtoseconds, cor- 
responding to the bulk E2 g optical phonon mode as dis- 
cussed below. It is to be noted that in addition to Zn 
and Cd samples, we tried to measure coherent optical 
phonons also in alkaline-earth metal Mg (0001), however 
we observed only very weak electronic response and no 
oscillations in the TR signal (not shown). The amplitude 
of the coherent optical phonon is significantly larger in 
Zn than that in Cd. This suggests that the interband 
electronic transition, which contributes to the imaginary 
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FIG. 2: Transient reflectivity change obtained for Zn and Cd 
at the lattice temperature 7 K. The inset shows the imaginary 
part of the dielectric function Im(e) for Zn (closed circles), Cd 
(open squares), and Mg (open triangles) (Ref. 27). 
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FIG. 3: Transient reflectivity change obtained in Zn at vari- 
ous lattice temperatures. The inset shows the Fourier trans- 
formed spectra of the time domain data. 



part of the dielectric function Im(e) as shown in the inset 
of Fig. 2, governs the excitation of the coherent optical 
phonon, although according to Eq. (3) the phonon am- 
plitude depends also on the deformation potential. In 
fact, at the laser energy of 1.55 eV, Im(e) shows a dom- 
inant peak for Zn, whereas it shows a weak shoulder for 
Cd. For Mg, however, Im(e) is almost zero, supporting 
the suggestion that this mechanism is indeed the main 
driving force. 

Hereafter we will focus on dynamics of the coherent 
phonon response and fast electronic response in Zn be- 
cause the signal amplitude is much larger than that in 
Cd and therefore it is possible to make precise analysis 
of the time domain signal. 



A. Coherent phonon response 

Figure 3 shows the time-resolved TR signal (AR/R) 
obtained on Zn at various temperatures at the constant 
pump fluence of 9.2/iJ/cm 2 . The frequency of the coher- 
ent optical phonon observed in the Fourier transformed 
(FT) spectra (see inset of Fig. 3) is 2.32 THz at 7 K, be- 
ing in excellent agreement with that of the bulk Ei g mode 
observed by Raman scattering. [^JH^ The Fourier trans- 
formed spectra exhibit a redshift of the peak frequency 
and broadening of the line width of the E^g mode as the 
temperature increases. The coherent phonon signal in 
time domain, not including ultrafast electronic response 
appearing at t < 600 fs, was fitted by a damped har- 
monic oscillator; Ae~ Tt cos(uJE 2g t + 4>o), an d in this way 



the amplitude A, the decay rate (damping) T, the fre- 
quency ujE 2g j and the initial phase 4>o are extracted. 



1. Temperature dependence of damping and frequency of 
the E2 g mode 

In general, damping of the coherent phonon is gov- 
erned by population decay (inelastic scattering) and/or 
pure dephasing (elastic scattering). In semimetal and 
semiconductor crystals, the decay process of the coher- 
ent phonon is dominated by the population decay due to 
anharmonic phonon- phonon coupling, rather than pure 
dephasing. [23, H(j The anharmonic decay rate depends 
strongly on the lattice temperature, while the pure de- 
phasing, e.g., e—ph scattering, depends in metals on the 
hot electron density. In order to examine the decay 
process of the coherent optical phonon, temperature de- 
pendence of the decay rate at the constant pump fluence, 
i.e., constant hot electron density, is measured. 

Figure 4 shows the decay rate of the coherent optical 
phonon as a function of the lattice temperature. The 
decay rate increases upon increasing the temperature. 
This behavior is well explained by the anharmonic de- 
cay model. |29| in which the optical phonon decays into 
the two acoustic phonons with half the frequ ency of the 
optical mode and with opposite wave- vectors, [29j 30] 
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FIG. 4: Decay rate of the Ei g mode together with the fre- 
quency (inset) as a function of the lattice temperature. The 
solid line represents the fit to the data using Eq. (4) . 

Here Tq is the effective anharmonicity as the fitting pa- 
rameter and ks the Boltzmann constant. Tq is deter- 
mined to be 0.06 ps _1 . The good agreement of the 
time domain data with the anharmonic decay model in- 
dicates that the damping of the coherent E 2g mode in Zn 
is due to anharmonic phonon-phonon coupling (popula- 
tion decay). The frequency of the E 2g mode decreases 
as the temperature increases as shown in the inset of 
Fig. 4. This temperature dependence is qualitatively 
in good agreement with the anharmonic frequency shift 
observed by Raman scattering measurements. |28| Such 
a frequency shift due to the lattice anharmonicity was 
also observed in III-V semiconductors, which was re- 
produced by ab initio calculations including various an- 
harmonic contributions (thermal expansion, third-, and 
fourth-order anharmonicity) . |3 1| 



2. Temperature dependence of the amplitude of the E^ a 
mode 

As shown in Fig. 5 (a), the amplitude of the coher- 
ent Ei g phonon significantly decreases upon increasing 
the lattice temperature. The decrease in the amplitude 
of the coherent optical phonon at higher temperatures 
is contrary to the temperature dependence of the Ra- 
man intensity. |2lL l32l | confirming that the amplitude of 
the noncquilibrium coherent phonon is not determined 
by Bose- Einstein statistics = l/[exp(hfl /kBT) — 1], 
which shows an increase in the phonon occupation num- 
ber at higher temperatures. The very close behavior of 
the temperature dependence of the phonon and electronic 
[see below in Fig. 9 (a)] responses suggests that the am- 
plitudes of the fast electronic and the coherent phonon 




FIG. 5: (a) Amplitude of the coherent Ei g phonon as a func- 
tion of lattice temperature, (b) Amplitude of the coherent 
E2 g phonon vs the pump fluences at the two typical temper- 
atures. The solid curves are the fit to the data using Eq. 
(12). 



responses have the same origin. We will discuss the tem- 
perature dependence of the amplitude of the E 2g mode 
in Sec. V. 



3. The photoexcitation dependence 

The photoexcitation intensity dependence of the co- 
herent E 2g mode was investigated in the range up to / 
= 18.4/iJ/cm 2 . The amplitude monotonically increases 
with the pump fluence both at 7 and 80 K as shown in 
Fig. 5 (b). On the other hand, the frequency and the 
damping of the coherent E 2g mode were independent on 
the pump fluence in this fluence range. The latter re- 
sult shows that the damping of coherent E 2g phonons 
in Zn is dominated by anharmonic phonon-phonon cou- 
pling, which is determined by the lattice temperature, 
rather than by pure dephasing via scattering by hot 
electrons. 16] 

When increasing the pump fluence by three orders of 
magnitude (up to 6 mJ/cm 2 ) using the amplified fem- 
tosecond laser system ( pulse duration of 130 fs) the 
time domain signal drastically changes. Figure 6 shows 
transient reflectivity change observed at 6 mJ/cm 2 . The 
width of the initial electronic transient response is dif- 
ferent from those observed in Figs. 2 and 3 because the 
pulse duration used in Fig. 6 is much longer than that 
in Figs. 2 and 3. Importantly, although the pump flu- 
ence increases by ~ 10 3 , the electronic response increases 
only by one order of magnitude, suggesting some kind of 
a saturation behavior. |33j The frequency of the coherent 
E2 g mode in the FT spectra (the inset of Fig. 6) is 2.10 
THz, which is slightly redshiftcd from that obtained at 
9.2/zJ/cm 2 at 295 K, 2.15 THz. This redshift would orig- 
inate from electronic softening of the latti ce l34| or ther- 
mal expansion due to lattice anharmonicity 35, 36] under 
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FIG. 6: Transient reflectivity change observed at 18 mJ/cm 2 
at 295 K. The inset shows FT spectrum of the time domain 
data, showing the coherent E^g mode. 
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FIG. 7: (Color online) Transient reflectivity change observed 
for Zn at 7, 100, and 295 K. The closed circles are experimen- 
tal data, and the solid curve is the fitting with Eq. (5). 



high-density photoexcitation. We note that strong back- 
ground signal from lattice heating arises at t > 2 ps and 
therefore the effect of thermal expansion may dominate 
the redshift of the coherent Ei g mode in Zn. 



B. Fast electronic response 

The subpicosecond electronic response measured at 
several temperatures is presented in Fig. 7. In order 
to fit the time domain data and to obtain the amplitude 
and relaxation time of the fast electronic transient, we 



utilize a linear combination of a damped harmonic oscil- 
lator and a single exponential decay function 



AR(t) 
R 



= H(t) [Aer rt cos(u) E2g t + <j) ) + Be^ t/Tq + const.} . 

(5) 

Here H{t) is the Heaviside function convoluted with 
Gaussian to account for the finite time-resolution, while 
B and r q are the amplitude and relaxation time of the fast 
electronic transient, respectively. In simple metals, such 
as Au, Ag,0 ultrafast transient is usually attributed to 
the electron-phonon thermalization, which is treated by 
the so- called two-temperature model (TTM) . The main 
idea of the TTM is that because of the fact that electron- 
electron thermalization is much faster than electron- 
phonon thermalization, electrons quickly (in tens of fem- 
toseconds) thermalize to a temperature T e , which can be 
much higher than the lattice temperature T; , resulting in 
Fermi-level smearing. In the next stage, electrons ther- 
malize with lattice in a characteristic electron-phonon 
thermalization time, which is in metals typically in the 
100 fs - 1 ps range. Since after tens of femtoseconds, 
when electrons have already thermalized, the changes in 
the occupied electronic density of states are limited to 
energies of ksTe near Fermi level, the photoinduced re- 
flectivity dynamics tracks the time-evolution of the elec- 
tronic temperature. The TTM is given by the set of two 
coupled heat equations |23, H3 



ciy = <?Cr e -r ; ), 



g(T e -T t ) + S(z,t), (6) 

(7) 



where C e (T e ) and C; are the respective heat capacities of 
electrons and lattice, g [= g{Ti)} is the e — ph coupling 
function and S(z, t) describes the absorbed energy, where 
z is the depth coordinate. In the pump-probe experiment 
S(z,t) has a Gaussian temporal shape. In the limit of 
weak perturbation (when change in the electronic tem- 
perature is small compared to the initial temperature) 
the relaxation of the electronic temperature is exponen- 
tial with e — ph thermalization time[j| 



-ph 



1 C e Q 
gC e + Ci 



(8) 



g(T{) is in the linear response limit particularly simple 
in the case of simple metals where the electron band- 
width is much larger than the Debye temperature &d, 
and the Debye model of the electron-phonon coupling can 
be used. In this case g(T) = dG(T)/dT, where |l|3a 



G{T)=A 9oo {-±-f 4^- 
6fl Jo e x - 1 



(9) 



Given that, the T-dependence of T e _ p /j Eq.jS] is com- 
pletely determined by the e—ph coupling constant (ffoo), 
Qd, C e (T), and Ci(T). Since the lattice specific heat 
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FIG. 8: (color online) Temperature dependence of the 
electron-phonon thermalization time (r q ). The closed squares 
are experimental data, whereas the solid curve is a fit to the 
data with the TTM model, described by Eq. (10). 



Ci is a factor 10 2 larger than the electronic specific heat 
(Ci » C e ) in a wide temperature range, and since the 
electron heat capacity C e = jT e ~ 7T; (when change in 
the electronic temperature is small compared to the ini- 
tial temperature), the e — ph thermalization time T e - p h 
is given by 

_ C e jJ\ 

T e-ph ~ trr\ ~ frn\i V iu J 

9{Ti) 9(Ti) 

For Ti < Q D the function g(T t ) varies Tf and for Ti > 
®d the function g{T{) becomes constant (goo)|2(- Since 
e D = 234 K for Zn, TTM predicts ~ T," 3 at T t < 

6_d/5 and r e _ p/l ~ 7} at T ( > 6d/5.|23 

The relaxation time of the subpicosecond transient in 
Zn obtained by fitting the data with Eq. (JjjJ is shown 
in Fig. 8 as the function of lattice temperature. The re- 
laxation time monotonically increases from 75±30 fs to 
260±30 fs as the temperature increases, similar to the 
behavior obtained for noble metals . The behavior fol- 
lows the prediction of the TTM at temperatures above 
50K, whereas no upturn in relaxation time is observed at 
low temperatures. The absence of upturn in relaxation 
time at low temperatures is common in metals, and is 
attributed to the fact that the basic assumption of the 
TTM that e — e thermalization time is fast compared 
to the e — ph thermalization fails at low temperatures. 
Therefore a nonthermal electron model (NEM) needs to 
be introduced that accounts for the low-temperature sat- 
uration of relaxation time [{J |39| - the analysis using a 
NEM is beyond the scope of this paper. By fitting the 
temperature dependence of the relaxation time with Eq. 
HI U| l we obtain the value of the e — ph coupling constant 
.9oc = 6.4 x 10 16 Wm _3 K _1 , which is the same order of 
magnitude as in other simple metals. |6j 
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FIG. 9: (a) Amplitude of the fast electronic transient as a 
function of temperature. Inset shows the corresponding ini- 
tial (T e - dashed line) and final (Ti - solid line) electronic 
temperatures as a function of temperature, (b) Amplitude of 
the fast electronic transient vs the pump fluences at the two 
typical temperatures. The solid curves are the fit to the data 
using Eq. (12). 



1. Temperature dependence of electronic transient 

As shown above, the temperature dependence of the 
relaxation time follows the prediction of the thcrmomod- 
ulation scenario, where the data above 50 K can be well 
fitted using the standard TTM. On the other hand, the 
amplitude of the fast electronic response, B in Eq. J^}, 
also shows pronounced temperature dependence as shown 
in Figure 9(a). Moreover the same temperature depen- 
dence is obtained also for the amplitude of the oscillatory 
transient due to the photoexcited coherent Ei g phonon, 
suggesting that the amplitudes of the fast electronic and 
the coherent phonon responses have the same origin. 

It should be noted that the temperature dependence of 
the amplitude of the photoinduced reflectivity transient 
has thus far not been studied for the case of simple met- 
als. Below we derive the temperature dependence of the 
transient amplitude B as would be expected in the case 
of the thermomodulation mechanism and show extremely 
good agreement with the data. 

To account for the temperature dependence of B in Eq. 
(5), we assume that the relaxation processes after excita- 
tion with short laser pulse follow the TTM, i.e., electron- 
electron thermalization is much faster than e — ph ther- 
malization. Furthermore, we assume that the changes 
in the reflectivity at optical frequency (1.55 eV) are due 
to photoinduced absorption, with photoinduced quasi- 
particles (in the energy range of fcsT e around the Fermi 
energy) as initial states and with final unoccupied states 
well above the Fermi energy. |4Qj, |4l| We should note that 
similar arguments have been successfully applied for the 
case of optical pump-probe spectroscopy in superconduc- 
tors and charge density waves. |43j Therefore, using 
the Fermi golden rule, the amplitude of the transient is 
proportional to the change in the occupied electronic den- 
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sity of states near Ep, i.e., the amplitude of the photoin- 
duced transient should be proportional to the photoin- 
duced quasiparticle density n p . In the limit when tem- 
perature is much lower than the Fermi energy, which is 
the case here, the number density of the electron-hole 
pairs n is in the Landau Fermi liquid exactly propor- 
tional to temperature, n oc T. We assume that after 
photoexcitation all the absorbed energy initially goes in 
the electronic subsystem, and the electrons are described 
by the increased electronic temperature T' e , Following 
this one obtains 

B oc n p — nr> — nr^ oc T e — T e , (11) 

where nx e is the quasiparticle density at the initial tem- 
perature T e , while tit' is the quasiparticle density af- 
ter electrons have thermalized to final temperature T' e . 
Taking into account that the electronic specific heat 
C e = 7T e , where 7 is the Sommerfeld constant and using 
the energy conservation law it follows 

B^T' e -T e = ^Tf+mh-T e , (12) 

where Ui is the absorbed energy intensity. 

There are several important implications of Eq. (|12fl . 
First, the amplitude of the transient is maximum at low 
temperatures and decreasing as the sample temperature 
is increasing. Second, at low temperatures the model pre- 
dicts a sublinear excitation intensity dependence, while 
at high temperatures the dependence becomes more and 
more linear. 

In Fig. 9 we show that the model, indeed, well repro- 
duces the observed temperature dependence. Here Ui/^y 
was taken as a fitting parameter since the estimated ab- 
solute value of Ui is subject to some uncertainty (possible 
errors in the spot size, reflectance, and the optical pene- 
tration depth may give rise to substantial error) . Impor- 
tantly the final temperature T' e obtained from the same 
fit [shown in the inset to Figure 9(a)] matches well the 
estimate given in Ref. 23. 

Although very good agreement with the temperature 
dependence of B is obtained at temperatures higher than 
~ 50 K, B shows saturation below 50 K, while the model 
predicts further increase in the amplitude. We ascribe 
this discrepancy to the fact that at low temperatures the 
TTM fails, since e — e thcrmalization time becomes com- 
parable or even longer than the e — ph thcrmalization 
time. [39]] 

Figure 9(b) presents the excitation intensity depen- 
dence of amplitude measured at 7 and 80 K, together with 
the fit using Eq. I|12|) using the same value of Ui/j as ex- 
tracted from the fit to the temperature dependence data. 
As seen, the model predicts the sublinear dependence at 
7 K, which is not observed in this intensity range, while 
good agreement with 80 K data is obtained. We believe 
that for the same reasons that the temperature depen- 
dences of relaxation time and amplitude at low temper- 
atures do not follow the TTM prediction; the fact that 
electrons are nonthermal on this timescale can give rise 



to the absence of the sublinear intensity dependence at 
low temperatures. 

V. DISCUSSION 

According to the model by Stevens et al. [2(j under 
the condition that the lattice temperature is constant, 
the amplitude of the coherent phonon depends on the 
imaginary part of the dielectric function Im(e) as well 
as the pump intensity I(t) [see Eq. (3)]. Here we re- 
view the temperature dependence of the coherent phonon 
amplitude in Zn in terms of Eq. (3). When the lattice 
temperature increases, thermal expansion occurs and the 
volume v c increases. This effect will increase the phonon 
amplitude according to Eq. (3). The phonon frequency 
decreases for sa 8 % with increasing the temperature from 
7 to 295 K, suggesting an amplitude increase [Eq. (3)] 
by wl6 % contrary to our experimental results, which 
shows a decrease in the coherent phonon amplitude as 
the temperature increases. The deformation potential 
So can be taken as nearly constant through the entire 
temperature range in the present study. [44| Thus, any 
parameters other than Im(e) in Eq. (3) cannot account 
for the temperature dependence of the coherent phonon 
amplitude. 

The initial phase of the coherent optical phonon is esti- 
mated to be 5° ± 5° at all temperatures from the fitting 
the oscillatory part with the single damped harmonic os- 
cillator, indicating that the coherent phonon oscillation 
follows cosine behavior. The cosine phase is consistent 
with the prediction from Stevens' model [see Eq. (2)], 
where the driving force of resonant ISRS has displacive 
character. Further argument for the displacive nature of 
the observed transient is the strong correlation between 
the photoinduced quasiparticle density and the ampli- 
tude of the coherent optical phonon observed in both 
temperature and fluence dependences. For example, in 
Fig. 5 we use the same fitting function as used for de- 
scribing the temperature dependence of the electric tran- 
sient. Fitting of the coherent phonon amplitude versus 
both the temperature and the fluence using Eq. (12) is 
surprisingly good. As discussed in Sec. Ill, the driving 
force F(t) oc J_ \E(t')\ 2 dt' is the Heaviside function 
convoluted with Gaussian and causes displacive behavior 
for the excitation of Raman active mode through reso- 
nant ISRS process. Because Im(e) is directly related to 
photo excitation of nonequilibrium electrons, quasiparti- 
cle density should follow Im(e): n p oc Im(e). Thus under 
the resonant condition, Qo oc Im(e) and n p oc Im(e). 

VI. CONCLUSIONS 

We have investigated ultrafast dynamics of coherent 
optical phonons and nonequilibrium electrons in transi- 
tion metal Zn and Cd using a femtosecond pump-probe 
technique. The optical response in time domain is well 
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fitted by linear combination of a damped harmonic os- 
cillation due to excitation of coherent E 2g phonon and 
a subpicosecond electronic transient. Dynamics of both 
the noncquilibrium electron distribution and the pho- 
toexcited coherent Ei g phonon were observed in a wide 
temperature range from 7 to 295 K. The amplitudes of 
both the coherent optical phonon and electronic transient 
show pronounced temperature dependencies. Based on 
the fact that the relaxation time of the fast electronic 
transient closely follows the prediction of the thermo- 
modulation scenario we analyze the relaxation dynamics 
in terms of the two-temperature model. Good agreement 
with the model is obtained between 50 K and room tem- 
perature, while below 50 K the model is found to fail 
similar to studies on simple metals, such as Au or Ag.[9J 
Below 50 K a nonthermal electron model |9j needs to be 
applied. 

In order to account for the temperature dependence of 
the amplitude we developed a model assuming the pho- 
toinduced absorption mechanism, where the signal am- 
plitude is proportional to the photoinduced change in the 
quasiparticle density. The model was found to account 
for the observed temperature dependence over wide tem- 
perature range. Importantly, the model predicts a sub- 
linear photoexcitation dependence at low temperatures 
which was not observed in the excitation range studied. 
This discrepancy may be due to the failure of the two- 
temperature model at low temperatures. The fact that 
the amplitude of the coherent Ei g phonon follows the 
temperature dependence of the amplitude of the fast elec- 
tronic transient suggests that under the resonant condi- 
tion both electronic and phononic responses are propor- 



tional to the change in the dielectric function. 

The damping and the frequency of the coherent E2 g 
phonon were also found to depend on the lattice tem- 
perature, which was explained by anharmonic phonon- 
phonon coupling, rather than pure dephasing, and by 
anharmonic frequency shift, respectively. There was no 
dependence of the frequency and the damping of the Ei g 
phonon on the pump fluence, suggesting that pure de- 
phasing due to scattering with hot electrons is negligibly 
small. When increasing the pump fluence up to the mil- 
lijoule per centimeter-squared range, the red-shift of the 
coherent Ei g mode frequency was observed, which was 
ascribed to thermal expansion. Since the Ei g mode in 
Zn may be sensitive to highly anisotropic compressional 
behavior of c/a axial ratio in Zn crystal (hep structure), 
as observed by Raman scattering under pressure, con- 
trolling the amplitude and the frequency of the coherent 
Ei g phonon will be useful for the study of the gradual 
anisotropic — > isotropic transition occuring in Zn.|45| 
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